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CHARGING CHARACTERISTICS OF MATERIALS: COMPARISON OF 

EXPERIMENTAL RESULTS WITH SIMPLE ANALYTICAL MODELS 

by Carolyn K. Purvis, N. John Stevens, and Jon C, Oglebay 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

An understanding of the behavior of materials, of dielectrics in 
particular, under charged particle bombardment is essential to the 
prediction and prevention of the adverse effects of sracecraft charging. 
This paper presents an effort to obtain such an understanding through a 
combined analytical and experimental approach. A one-dimensional model 
for charging of samples in the LeRC test facility is used in conjunction 
with experimental data taken in this facility to develop "material 
charging characteristics" for silvered Teflon. These characteristics 
are then used in a one-dimensional model for charging in space to ex- 
amine expected response. Relative charging rates as well as relative 
charging levels for silvered Teflon and metal are discussed. 


1. INTRODUCTION 

Two previous papers^ »2 described the test facility, test methods 
and measurements, and the results of various materials characteristics 
tests performed at the Lewis Research Center in support of the space- 
craft charging investigation. The present paper summarizes the analyt- 
ical work which has been performed interactively with this experimental \ 

work. The goal of the analysis is twofold. First, is to model the ( 

charging of material samples in terms of the material's parameters. 

Second, since a goal of the entire study is to predict behavior of 
spacecraft surfaces, an attempt is made to "scale" the environment, 
that is, to relate results obtained using a monoenergetic beam in the 
ground test facility to expected results with distributed particle 
fluxes of the space environment. An approach to establish this envi- 
ronment scaling is to develop models of charging for both charged par- 
ticle environments, and assume that the material properties are con- 
stant. Then differences between material charging behavior underground 
test and in space are a result of the differences in the two environ- 
ments. It is recognized that the vacuum levels in the ground test and 
space environments are also different. No attempt is made here to 
account for this p actor. 

^Published in the proceedings of the Conference on Spacecraft Charging 
Technology, NASA TM X-73537, 1977. 
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This paper, then, represents a first attempt at attaining the two 
goal9 of characterizing material charging and scaling to the space en- 
vironment, The models used are one-dimensional, and describe charging 
of samples in terms of the charging of a capacitor. This type of model 
has been used by a number of workers3«4>5,6 to describe spacecraft 
charging. 

The procedure used herein was to first develop a one-dime '.slonal 
model to describe charging of samples in the LeRC vacuum test facility. 
This model contained a number of parameters which were varied to pro- 
vide best fits to experimental data obtained in the facility. The val- 
ues of these parameters which yielded the best fit were Identified as 
the "material charging characteristics." These were then used in con- 
junction with a one-dimensional model for charging in the space sub- 
storm environment to make some predictions of the charging behavior of 
the materials in space. The insulator studied here is 5 mil silvered 
FEP Teflon kS/ 


2. ONE-DIMENSIONAL (1-D) GROUND TEST FACILITY MODEL 

In the LeRC test facility, a monoenergetlc beam of electrons with 
energies of to 20 keV is directed at normal incidence to planar sam- 
ples. A beam current density of 1 nA/cm2 was used to obtain all test 
data discussed in this report. 

The Ground Test Facility Model is a quasistatic current balance 
model. The current densities considered are those due to primary (beam) 
electrons, secondary electrons, backscattered electrons, and leakage 
current through the sample bulk. These are denoted by j e , j s , jgg, 
j respectively. The sample is assumed to charge like a capacitor. 
Tnus, a time balance equation is of the form 


dv s 

c dt + 1 z 


3 e “ ^s ^BS 


CD 


where Vg is the magnitude of the surface voltage. All signs are ex- 
plicit in this equation, and in all others used in the test facility 
model, that is, all symbols stand for positive numbers. The actual 
surface voltage is of course negative, and this is assumed throughout. 
Thus, this model does not predict positive surface voltages correctly, 
since positive surface voltage would tend to reattract the emitted 
secondary electrons and this effect has not been included here. 

The first term in equation (1) represents the net charge deposited 
on the surface (per cm^) in a time step, and so is termed the charging 
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current density, and denoted by j c . The procedure used was to ini- 
tialize at t « 0 with Vg “ 0. Then AVg was calculated from 


AV i r -i 

1t Sa U j e- Js “ JBS -JJ < 2 > 

The current densities are all functions of Vg. Their functional forms 
are given in figure 1. Equation (2) is solved by an iterative proce- 
dure. This equation can be expressed as 


AV s = f (V g) At 


(2a) 


With the initial assumption Vg = t = 0, a suitable At is chosen, and 
AVg calculated. Vg is then incremented by AVg (set = AVg for the 
first iteration) ana the procedure repeated until equilibrium is reached, 
that is, until - C dVg/dt = 0. 

The several current densities in the above equations are functions 
of surface voltage. All but the leakage current density are functions 
of the primary electron beam voltage and current density. The param- 
eters which can be varied are the secondary emission maximum yield, 6 m , 
and energy for maximum yield, V , the backscatter coefficient, £, the 
resistivity, p, and the capacitance C. In practice, values for S ra 
and V m were taken from the literature; values for p were determined 
from the measured surface voltages and currents at equilibrium, and fits 
obtained by varying % and C. This is discussed more fully in sec- 
tion 4 below. 

It should be noted that this model does not account for beam spread 
in angle or energy, the presence of the vacuum tank walls, or residual 
gas in the chamber. Since it is one-dimensional, it cannot, of course, 
explicitly describe edge effects, or effects due to surface variations. 

Derivations of the current density equations are presented in 
Appendix A; a summary of the model and equations is given in figure 1. 


3. COMPARISON WITH EXPERIMENTAL DATA 
3.1 "Floating" Aluminum Plate 

The first test data used in conjunction with the test facility 
charging model were those taken with a bare aluminum plate which is nor- 
mally used for substrates. This plate was mounted in front of a second 
identical pL.ce which was grounded to the chamber walls. These two 
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plates were held apart and electrically isolated by a Teflon spacer plug 
0.7 cm long. The plates are rectangular, with dimensions 15 cm x 20 cm 
(-300 cm^ area). Surface voltage of the floating plate was measured as 
a function of time using the surface voltage probel. 

Figure 2 shows these data and the best fit calculated curves. In 
obtaining these fits to the data values of 6 m = 3 and Vm a 400 were 
used. These are consistent with ranges of values for these parameters 
given by Gibbons^ for a surface layer of AI 2 O 3 . It is reasonable to ex- 
pect some oxide on the surface of the "bare aluminum" plate since it had 
in fact been exposed to air. This points up the strong dependence of 
charging phenomena on sample surface condition, and indicates that care 
must be taken in making predictions for charging to consider the state of 
the surface. 

With these values in the expression for secondary emission, the val- 
ues of C used to obtain the curves in figure 2 varied from 1.2xl0 _ ^0 
farads at Vg = 5kV to 5x10“ ^ farads at Vg = 16 kV, decreasing approx- 
imately linearly with Increasing beam voltage. The values used for the 
"backscatter coefficient" varied from 0.15 at Vg - 5kV to 0.5 at 
Vg = 16 kV, again in approximately linear fashion. The expression used 
to calculate secondary emission current density is derived from an ex- 
pression for yield as a function of primary energy due to Sternglass^. 

It is plausible that the required variation of £ to obtain fits to the 
data is accounted for, in part, by deviations of the actual secondary 
yield from that predicted by Sternglass formula. That is, the adjustments 
to X, represent adjustments to the sum of the backscattered and secondary 
electrons. 

The final point of concern here is the time scale for charging of the 
floating plate; it reaches its equilibrium floating voltage with a time 
constant of several seconds. This is not surprising, since the capaci- 
tance of the test plate to its surroundings is expected to be small... The 
time scale is relevant, however, to the question of the behavior of com- 
posite samples. This is discussed more fully In section 3.3 below. 


3.2 Silvered Teflon Samples 

The 1-D model was next used to fit current and surface voltage data 
from silvered Teflon samples. These consist of an aluminum substrate 
with three strips of 5 mil silvered Teflon mounted with conductive adhe- 
sive to the substrate. Each of these strips was 5 cm wide and 20 cm long. 
During test, the aluminum substrate (and consequently the silver) was 
grounded, while the Teflon surface was bombarded with electrons. 

The data and calculated fits for beam voltages at which equilibrium 
is reached are shown in figure 3. These data are a composite of four 
separate data sets, and indicate charging times on the order of minutes. 
The error bars reflect the scatter in the data as well as the +5% 
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resolution uncertainty in the voltage measurements. Since for insula- 
tors there are strong voltage gradients near the edges of the samples, 
the surface voltage measurements are those read at the sample center 
which is uniform. 

To obtain these curves, the effective resistance was calculated 
from the equilibrium values of surface voltage and leakage current. 

These values indicate an effective resistivity for the samples of about 
9x10^5 fi-cm, about an order of magnitude below published values for 
Teflon. If one assumed that the published value of lO-^ fi-cm is accu- 
rate, then a parallel path having a resistance of about 4xl0^ fl is in- 
dicated by the data. This could be a surface leakage, or an edge leakage, 
or leakage through the sheath. 

Values of <S m = 3 and V m = 300 were used in the expression for 
secondary emission, in accordance with the data given by Willis and 
Skinner^. The values of C required ranged from 14.6 pf/cin^ at 
Vg «* 5 kV to 10 pf/cm at Vg = 12 kV and decreased monotonically 

with increasing beam voltage. Extrapolation of the curve to Vg = 0 in- 
dicates a dielectric constant for Teflon of '2.1. This decrease in ef- 
fective capacitance is believed to be associated with edge effects. Edge 
gradients are observed in the data. These become more pronounced at 
higher voltages, reducing the effective area^. The value used for the 
backseat ter coefficient varies about an order of magnitude for the Teflon 
samples. Not only does it change with beam voltage, but also appears to 
change during charging at a single beam voltage. For the initial portion 
of the charging curves, C varied from 0.25 at the lower beam voltages to 
0.05 at high beam voltages. At equilibrium, C varied from 0.2 at low 
beam voltages to 0.25 at high ones. These variations were not linear; 
rather, C was relatively constant at low and high beam voltages, with a 
transition occurring between V B *= 8 kV and Vg = 12 kV. Again, part of 
these variations may be due to deviations of the actual secondary emis- \ 

sion from that calculated by the analytical expression being used. \ 

Some investigation was undertaken to study the behavior of the 
Teflon samples during arcing. Figure 4(a) shows a curve fit to a com- 
posite of two data sets for the initial charging transient with a beam 
voltage of 16 kV. Figure 4(b) shows the same calculated curves, this 
time with a single set of data, The left hand set of curves and points 
are the initial charging transient. An arc occurred on this sample be- 
tween the time of the voltage reading at t = 4 min and that of the cur- 
rent reading taken at t - 5h min; the surface discharged, and charging 
was repeated. The curves shown for the post-arc charging transient are 
identical to those for the Initial transient, but shifted in time. This 
indicates that the charging transient is quite repeatable, not only from 
test to test of Teflon samples, but also through arcing. That is, at 
least short term, the arcing does not affect the charging behavior of 
Teflon samples. 



6 


h 


3.3 Composite Sample 

One of the stated objectives of the present work was to Investigate 
the "environment scaling" effects between the ground test and space en- 
vironments. An obvious difference other than the environment between 
ground tests so far described and the space condition is that the 
studies of silvered Teflon discussed above were all conducted with the 
substrate grounded. In contrast, for the case of a spacecraft in orbit, 
the entire body, including the "grounds" must come into equilibrium with 
the charged particle environment. The question of the behavior of a com- 
posite system becomes particularly interesting because of the divergent 
time scales in which charging of "floating" metal plates and silvered 
Teflon above ground are observed to occur. 

To investigate this question, a composite sample was built and 
tested. The sample consisted of an electrically floating standard alum- 
inum substrate with two strips of the 5 cm silvered tape mounted on it. 

This allowed for a 5 cm strip of the aluminum between the two Teflon 
strips to be exposed to the electron beam. The aluminum substrate was 
mounted in the test chamber in the same manner as the floating aluminum 
plate described in section 3.1. This configuration is depicted at the 
top of figure 5(a). The bottom of this figure shows a voltage trace at 
equilibrium for thiB sample. 

The expectation was that this composite sample would charge in two 
stages because of the different effective capacitances through which the 
aluminum and the Teflon must charge. The prediction, shown in figure 5(b), 
is based on the idea that when the beam is turned on, the aluminum should 
charge to its equilibrium voltage with its time constant of seconds, car- 
rying the Teflon voltage with it. When this has occurred, the Teflon 
should continue to charge from the equilibrium voltage of the aluminum to 
its own equilibrium voltage with its own time constant, that Is, minutes. 
The curves in figure 5(b) were thus obtained by superimposing the curves 
for aluminum alone and for Teflon alone with the Teflon curve shifted, so 
that it coincides with the aluminum curve at the point of equilibration 
for the aluminum. As can be seen from the data plotted in this figure, 
the expected behavior was found. 

The expectation that the aluminum and Teflon comprising the compos- 
ite sample should charge to the same surface voltages as had the float- 
ing plate and the Teflon samples above grounded substrates was based on 
the observation that in both those cases the equilibration was dominated 
by secondary emission phenomena rather than by leakage currents to 
ground. This is evidenced in two ways. First, plots of surface voltage 
at equilibrium versus beam voltage are straight lines; such behavior is 
supposed to be associated with emission dominated equilibration. Second, 
examination of printouts of the model calculations reveals that, at equi- 
librium, the leakage current density term is several orders of magnitude 
smaller than the other currents in the model. The conclusion, then, is 
that for this type of composite sample, each part responds to the 
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charging environment with its characteristic time constant, and comes 
into equilibrium at its characteristic surface voltage so long as leak- 
age current does not play a dominant role in the equilibration. 

It should be noted here that the tests run on this composite sample 
were not extensive. Further experimental investigation of this and other 
composite samples are planned. 


4. ONE-DIMENSIONAL SPACE CHARGING MODEL 

In order to predict charging behavior in space, a one-dimensional 
model for charging in space analogous to th* one-dimensional ground test 
model was developed. The essential dlfferei.ee between the two models, 
is that, the space model assumes an Isotropic Maxwellian particle dis- 
tribution containing both electrons and ions, and a spherical collection 
geometry. The current densities are derived from Langmuir probe calcu- 
lations. This type of calculation has been used by several authors^ ,5, 6, 10 
to treat the spacecraft charging problem. Derivation of the current den- 
sity equations is given in Appendix B; a summary description is presented 
as figure 6. As 1b indicated in this figure, the model as presented and 
used here assumes a geomagnetic substorm condition. That is, it assumes 
that Vs is negative, so that electrons are repelled and ions attracted. 

In these equations Vg is an algebraic quantity, that is, the sign is 
Implicit, 

As is evident from the sketch in figure 6, secondary electrons due 
to ion impact are not accounted for in this model. This is because the 
Intent here is to use material charging characteristics found by fitting 
the ground test model to experimental data In conjunction with the space 
mode] to predict space charging. Since there were no ions used in the 
experiments, no coefficient for secondary electrons due to Ion impact was 
determined. Therefore, this current density source is not considered. 

The surface voltage values predicted by this model are therefore somewhat 
larger than if secondary electrons due to ion Impact had been included. 

For example, a secondary coefficient of 1 would yield about 10% reduction 
in the equilibrium voltage calculated for aluminum at V e - 5 kV. 

The procedure used to calculate charging is identical to that de- 
scribed for the ground test model. Material characteristics used were 
those determined by fitting the ground test model to the data. The 
capacitance of the aluminum (considered to represent the spacecraft 
"ground") was taken to be 15xl0”H farads, which is the capacitance of a 
1 meter diameter sphere (to represent a "typical" spacecraft dimension) 
to infinity. This capacitance was chosen because the relevant capaci- 
tance for charging floating metal objects is that of the object to its 
surroundings. Those parameters (notably C and £ for Teflon, and £ 
for aluminum) which varied as functions of beam (and therefore surface) 
voltage were associated with the equilibrium surface voltage for the 
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appropriate test for purposes of making the space voltage calculations. 
The relationships between electron and ion temperatures and between 
temperatures and current densities were taken from the Provisional 
Specification for the Geomagnetic SubBtorm Environment. 11 This Speci- 
fication is given as figure 7. Thus, results of the space calculations, 
shown in figure 8, ore given as functions of electron temperature only. 

In figure 8(b) two curves are shown for the surface voltage of sil- 
vered Teflon as a function of electron temperature. The first curve 
calculated used the experimentally determined value of 9x10^ fi-cm for 
the effective resistivity of Teflon. This curve bends sharply to the 
right as electron temperature increases. An inspection of the current 
densities driving the equilibrium indicated that leakage current played 
a large part in the equilibration of the Teflon. Thus, this curve 
yields a "good" value for the -urfaco voltage of Teflon if the space- 
craft "ground" is actually ne„ plasma ground. This would be the case 
if, for example, photoemission were holding the spacecraft ground near 
plasma ground and the Teflon surface of concern were shaded. However, 
if the spacecraft is assumed to be in eclipse, the aluminum representing 
the spacecraft ground is predicted to acquire a large negative voltage 
(solid line). In this case, it is clear that leakage current cannot 
drive the equilibration of the Teflon surface according to the dashed 
line. Therefore, a second curve was calculated for the Teflon, based 
on the assumption that there was no leakage. This is shown in the dash- 
dotted line. 

Figure 8(a) shows the rates of charging predicted for the space 
condition. The time to charge the aluminum is predicted to be seconds, 
comparable to the time required to charge the floating plate in the 
ground experiments. The time required to charge the Teflon is predic- 
ted to be several minutes at lower electron temperatures, ranging to 
tens of minutes as the electron temperature increases. Thus, equilibra- 
tion of Teflon is predicted to require significantly longer time in 
space than it does in ground experiments. 

Using the curves shown in figure 8, it is possible to predict the 
response of a "spacecraft," composed of an aluminum structure partially 
covered with silvered Teflon, to substorm and eclipse conditions. Such 
a set of predictions is shown in figure 9. For purposes of this figure 
it is assumed that photoemission is sufficient to hold illuminated sur- 
faces close to plasma ground. 

Figure 9(a) assumes that the Teflon surface of interest is shaded 
when the spacecraft is in sunlight. A substorm injection with a 5 keV 
Maxwellian electron distribution is assumed to occur at t = 0 and this 
environment is assumed to remain constant throughout the time shown. An 
eclipse is assumed to occur from t = 60 minutes to t = 120 minutes. 
When the substorm occurs, the dark Teflon charges according to the 
dashed curve of figure 8(a), shown in figure 9 as a solid curve; the 
aluminum at "ground" is assumed to be held near plasma ground by 



photo emission. The Teflon surface reaches Its leakage dominated equi- 
librium voltage of -8.5 kV with Its rime constant of about 20 minutes. 
When the spacecraft enters eclipse, the aluminum charges quickly (in 
seconds) to its equilibrium voltage of -10 kV. It is assumed that the 
charge on the Teflon surface is immobile on this time scale, so this 
surface remains at -8.5 kV. At this point, the leakage current, which 
was driven by the voltage differential of 8.5 kV disappears, so the 
Teflon finds Itself to be no longer in equilibrium with its environment 
and proceeds to charge to its "floating” value of “-14.5 kV in a char- 
acteristic 20-30 minute period. Upon exit from eclipse, a similar pat- 
tern is followed. The aluminum falls quickly to near plasma ground. 
Because this discharging is driven by photoemission, it requires only 
about 0.02 second for the aluminum to reach plasma ground (assuming 
-10"9 amps/cm 2 photocurrent). Now, the Teflon again finds itself out 
of equilibrium with its plasma environment, and proceeds to discharge 
slowly to its previous equilibrium potential of -8.5 kV. 

Figure 9(b) shows a similar type of time history for an Insulating 
surface which is exposed to sunlight. Again, the solid line represents 
the surface voltage of the Teflon and the dashed line the spacecraft 
ground. The entry into eclipse and subsequent charging up is analogous 
to the charging of the composite sample discussed in section 4 above. 

The aluminum charges rapidly (in seconds) to its equilibrium value. Be- 
cause the Teflon had no significant charge on its surface, its voltage 
follows that of the aluminum until the aluminum reaches equilibrium. 

The Teflon then continues to charge slowly to its equilibrium potential. 
Upon exit from eclipse, both the aluminum and the Teflon are discharged 
by photoemission. Thus, the aluminum reaches a plasma ground in about 
0.02 second, as in the previous case. The Teflon also discharges more 
quickly than it charged; it requires about 4 minutes to reach plasma 
ground . 

These results indicate a need for charging studies which take into 
account relative charging rates as well as different equilibrium charg- 
ing levels of various spacecraft surfaces. A "typical" spacecraft has 
several different types of surfaces (solar cells, thermal blankets, etc.) 
each of which can be expected to charge with its own time constant. The 
Importance of the effect of the different time constants should be 
assessed. 


5 . CONCLUDING REMARKS 

The present study has resulted in the development of a set of "mate- 
rial charging characteristics" which describe the charging of small 
(300 cm 2 ) samples of 5 mil silvered Teflon and oxidized aluminum. Based 
on these characteristics, predictions of charging in space have been made 
and used to estimate the behavior of a composite body under conditions of 
substorm and eclipse. Several Interesting differences between charging 


behavior under ground test conditions and predicted behavior in space 
have been noted for silvereO Teflon. Under ground test conditions 
silvered Teflon acquires the same surface potential mounted on grounded 
and floating substrates because the equilibration is douxnated by sur- 
face emission phenomena. However, in space, significant differences 
are predicted in surface voltage for these two mounting configurations. 
This is because leakage current dominates the equilibration of the 
Teflon when the aluminum is at ground, while surface emission dominates 
for the floating aluminum case. The time required for silvered Teflon 
to charge to equilibrium in tests is several minutes; this time is pre- 
dicted to be several tens of minutes in space. 

In contrast to the several minutes to several tens of minutes time 
scales for silvered Teflon, floating aluminum samples are observed in 
ground test and predicted in space to charge to equilibrium in seconds. 
Their capacitance is much lower than that of Teflon since it is determ- 
ined by their surroundings. There are also orders of magnitude differ- 
ences in time scales for discharging by photoemission of aluminum and 
Teflon. This discrepancy in charging and discharging rates gives rise 
to sudden changes in the electric fields which the Teflon must sustain 
upon entry into and exit from eclipse. It is felt that these differen- 
tial charging rates as well as differential charging levels may be im- 
portant and should be investigated further. Thus, transient (quasi- 
static) as well as steady state models should be developed for charging. 

The one-dimensional models described herein have been found useful 
in the interpretation of experimental results, and as guides to relating 
test results to expected space behavior. Models of ground test situa- 
tions are needed since they can be used interactively with test data. 
This is especially true since it is impossible, or at least impractical, 
to simulate accurately the geosynchronous environment. Thus, environ- 
ment scaling must be done through use of models, at least for the 
present . 

Finally, higher dimensional models are needed. One-dimensional 
models cannot account for such things as edge effects or interactions 
between adjacent surfaces at different potentials as with different 
charging properties. Such effects are clearly important, 2 and may 
dominate the charging behavior of multi-surface samples and spacecraft. 
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APPENDIX A 

ONE-DIMENSIONAL GROUND TEST MODEL 
A.l. INTRODUCTION 

In this model, electrons from the electron gun are assumed to ap- 
proach the planar sample normally. All motion is restricted to the 
x direction (see fig. A.l). The electron beam Is assumed momoener- 
getlc, with energy 

A tr 
e b * eV B 


where e is the electronic charge and V B the beam voltage. The cur- 
rent density emitted from the gun is given by 



where N 0 B particle density. It is assumed here that N 0 e (the charge 
density) is constant. Thus, the continuity equation requires that some 
particles are '‘lost." 

We wish to calculate current densities to the sample surface. Cur- 
rent densities to be considered are those due to primary electrons, sec- 
ondary electrons, backscattered electrons and leakage through the bulk 
of the Insulator. The Insulator is assumed to be mounted above a 
grounded substrate for purposes of calculating leakage. 

Throughout this development the sample surface is assumed negative, 
and all secondary and backscattered electrons are at.sr*"rd to escape. 

All signs are given explicitly, so that, symbols represent positive 
quantities. 


A. 2 CURRENT DENSITY DUE TO PRIMARY ELECTRONS 

Consider electrons approaching the sample with energy e . It is as- 
sumed that the electrons do not collide with one another, and that motion 
is restricted to one dimension. If the sample surface has a repulsive 
potential of magnitude Vg, energy conservation requires that they arrive 
at the surface with energy e-eVs» Since all electrons leave the gun 
with energy eV 0 , the current density to the sample surface is simply 
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Rearranging and using equation (Al) , we find 

\l/2 




/ Vg r 


(A3) 


A. 3 CURRENT DENSITY DUE TO SECONDARY ELECTRONS 

Sternglass 8 has given the following expression for secondary yield 
as a function of primary electron energy an impact; 
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where <$ m is the maximum yield, eV is the primary energy for which 
maximum yield is attained and is primary energy at impact. From 
the discussion given In (A2) above, for this case 


Si = eV B “ eV S 


<A5) 


The secondary current density is then 
j 8 = j e 6(«V B - eV s ) 


<A6) 


which is 
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A. 4 CURRENT DENS.ITY DUE TO BACKSCATTERED ELECTRONS 

No analytical expression was found for backscattered electron 
emission. For simplicity} it was therefore assumed that backscattered 
electron current density represents a fraction of the incident current 
density. Thus, 


^BS ~ C ^e 
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A. 5 LEAKAGE CURRENT DENSITY 
Leakage current Is generally represented by 

V S 

ia B X 
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In terms of bulk resistivity this is 


VsA 
H = p Z 


(A10) 


where p is bulk resistivity, A is the area and % the thickness of 
the sample. Then 


h 



(All) 


A. 6 THE l-D MODEL 

The primary electron current density represents a source of elec- 
trons arriving at the sample. The other three current densities repre- 
sent loss of electrons from the surface. Thus, the net current density 
to the surface is 


“ j e " j s " ^BS ” ^ A ^ AL 

iHs net current density plays the role of a charging current to the 
surface. Thus, if we represent the sample's charging as the charging 
of a capacitor, we have 


dVc 

= c “dt~ = je “ jg “ JbS ■ H ( A13 ^ 

where C is capacitance, here expressed in farads per square centi- 
meter to maintain consistency of units. Equation (A13) is solved in 
the manner described in the text on a computer to calculate the 
charging. 

It remains to associate the experimentally measured parameters, 
surface voltage, and total current to ground, with calculated values. 
The surface voltage association is trivial; it is simply -Vg. The 



total current to ground is the charging current plus the leakage cur- 
rent for the sample as a whole. Thus* 
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APPENDIX B 

ONE-DIMENSIONAL SPACE SUBSTORM MODEL 
B.l INTRODUCTION 

The one-dimensional space model assumes a two-dimensional iso- 
tropic Maxwellian velocity distribution for primary particles, and a 
spherical collection geometry. The calculations are essentially those 
for a spherical Langmuir probe. The present calculations are based on 
the work of Langmuir^ and Grard, at al.13 and follow closely the der- 
ivation of Cauffman.13 The latter work has not been published; there- 
fore, portions of it are reproduced here for clarity. Such portions 
are identified by superscript reference. 

10 

Geometry for the calculations is depicted in figure B.1. The 
sheath is assumed to have radius a, the collector radius R. The 
radial and tangential velocity components in the "undisturbed” region 
(sheath edge) are v r and v^., respectively and those at the collector 
surface u r and u t . The surface potential of the collector is Vg. 
The potential in the sheath is assumed to be a function of radial dis- 
tance from the collector and to be monotonic. The plasma is assumed 
collisionless, that is, orbit limited theory applies, and energy and 
angular momentum are assumed constant for each particle. 

The integral requiring solution for current densities due to pri- 
mary electrons^ >13 and ions, and baekscattered electrons is 
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where the subscript i stands for either e (electrons) or p (H+ ions), 
and 0r and er are the angle and energy at impact on the collector 
(see fig. Bl). 

For secondary electron current, solution must be found for 13, 13^ 

ir/2 

d e 4ira3 / d ^ e 

J s " 4irR 2 ~ 4itR 2 de R / d ^R de R (eRjBR)'S(eR) (B2) 

4 -'0 Uo 


where <5 (e r ) is the secondary electron yield as a function of electron 
impact energy. 
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These integrations cannot be performed directly because the dis- 
tribution function for the particles at the collector is unknown, and 
therefore we can not determine dj^/de^. However, we do know the dis- 
tribution at the sheath edge, and can therefore determine dji/de a . 

If we assume the plasma to be collisionless, we can also convert the 
limits on er and 0r to limits on e a and 0 a > and perform the re- 
quired Integrations on these variables. 


B.1.1 Conditions for Collection 

In order to contribute to current collected at R, a particle 
must have energy e R _> 0 and direction 0 £ 0^< ir/2. Since the 
plasma is assumed collisionless and V(r) is assumed monotoMc, each 
particle's energy and angular momentum must be conserved. Assume the 
particles of interest have charge -e. Energy conservation demands 

j Mi(u| + ut) ” \ MiCvr + Vt) + eVs <=^ e r “ e a + eV s (B3) 


Angular momentum conservation demands 


Ru t = avt^ Ru sin 0r = av sin 0 £ 


(B4) 


where 


u = (u| + u?) 1/2 and v = (vr + v£)^ 2 


Solving for e a and 0 a in terms of e^ and 0 R yields the condi- 
tions for collection 10 


0 < 0„ < £ => 0 < 0. < sin -1 !- 
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> e R > 0 => « > t a > e 0 = 


0 for Vg _> 0 (attraction) 


-eVg for Vg <_ 0 (repulsion) (B5) 
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B.1.2. Energy Format 


An Isotropic Maxwellian velocity distribution in three dimensions 
has a distribution given by 
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We are interested in a two-dimensional distribution which can be found 
from 


2ir 


g(v r ,vt) = f (v r ,v t sin ip t v t cos i|t)di|/ 



(B7) 


Substituting for f(v) and integrating, we have 


v3/2 
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Langmuir^ gives the incremental current across the sheath as 
di ° 4iTa^Nev r g(v r ,vt)dv r dvt 

Then, 

dii a 2 Nte / Mi \ 3/2 f M^vJ + v^)] 

d ^i = 4^ = i 2 727 Spi) VrVt eXP L 2eVi J dv * dv t 


(B9) 


(BIO) 


and, changing to e a , 0 a coordinates we have 


10 


a/2 

e 

\ 2 .... 

a 

5 e Q exp 
/ (eV.) 2 3 

“ eV ± 


sin 0 Q cos 0 & d0 a de a 


as required. 
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B.2 PRIMARY AND BACKSCATTERED PARTICLE CURRENT DENSITIES 

Since backscattered electron current density is considered to be 
simply a fraction of the incident electron current density, the same 
integration applies to both. The calculation for positive ion collec- 
tion is the same as for electron collection with approprlte sign changes 
to account for the positive charge, and using the ion mass and tempera- 
ture. In equations (B9) to (Bll), the sign on the leading charge has 
been suppressed. Appropriate signs will be supplied in section B.5 
below. 

The integration to be performed is 



(B12) 


Cauffman 1 ^ 


evaluates this 


integral and finds, 


1 *) 

in agreement with Langmuir . 
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Since the interest here is in modelling charging in substorms with no 
photoemission, we expect V s negative. Thus, electrons are repelled 
and ions attracted. So we have for electrons 



(B14) 
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for ions 


J 


P 


^PoV 1 ' 
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where the minus sign reflects the fact that ionB are attracted by nega- 
tive Vg. For backscattered electrons, 


hs 55 « e 


5Jeo 
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where ? is the backscatter coefficient. 


B. 3 SECONDARY ELECTRONS DUE TO ELECTRON IMPACT 

Secondary yield as a function of electron Impact energy has been 
given by Sternglass 8 as 


«(er) » 7.4 6 m exp 




(A4) 


where 6 m is the maximum yield and eV m the energy at which the maxi- 
mum yield is obtained. To determine secondary electron current density, 
we must multiply the left hand side of equation (Bll) by 6 (er) and in- 
tegrate. Thus, we need 
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which is Cauffman's^ 8 equation for secondary electrons, except that he 
uses a sec G r dependence of <S on 0 r which is not used here. The 
0 a integral is the same as before and yields 
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Now, sequentially setting x^ = e Q + eVg and n “ (eV e )“^ 2 x + (V e /V m )^^ 


and substituting, we find 


V e /V S V e 
J s = 2 J eo 7 * 4 6 m V m exp W a + V m 


exp(-n 2 )dn 
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Since we are considering substorm cases only here, the condition V g <_ 0 
is of interest. For this case the integral in (B19) is just 


V^j I exp(-rT)dn 


(B20) 
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Recalling the expression for repeated integrals of the error function 
complement-*-*'' 


i n erfc(z) - 7* 


[t - z 


n. r exp(-t 2 )dt 


(B21) 


where, by definition, 


i n erfc(z) 


L n “^erfc(t)dt 


i°erfc(z) = 7jf 


exp(-t 2 )dt = erfc(z) 


and identifying t, z, and n with the appropriate variables in (B20), 
we have 
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So, from (B19) and (B22) we have 


U “ G W-4 6 m v?5 l i 5 erfc 


fc) f v s v e 
to 6XP W + 


(B23) 


and we note that the dependence of j s on Vg is the same as that of 
j e and j BS . 

B.4 LEAKAGE CURRENT DENSITY 


Leakage current density is defined in the same manner for the space 
model as it was for the test facility model (see section A. 5 above). 
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Thus, we have 

, vs 

•H " p z 


(B24) 


where p Is bulk resistivity and l Is the thickness of the insulating 
film. 


B.5 THE 1-D SPACE MODEL 

Now, the net current density to the sample surface Is 

Jc " “Je + Jp + + ^BS " h (B25 > 

where the signs on the current densities are given explicity here, and 
we recall that Vo in this model is algebraic (l.e., can be positive 
or negative), although the derivations have assumed it negative. 

The net current density plays the role of a charging current to 
the surface, so that 

dVg 

j c » c -ar - -j e + Jp + + hs - h < B26 > 


where C is capacity expressed in farads per square centimeter and 
where we assume we are charging a capacitor. Equation (B26) is solved 
in the same way as equation (A13) on a computer to determine Vg versus 
time for charging. 


CURRENT DENSITIES TO SAMPLE 


£ ASSUMPTIONS 

^ MONOINERCETIC DIRECTED BEAM 

^'-CAPACITANCE -C NORMAL INCIDENCE ON PLANAR SAM PUS 

^SURFACE AT 
VOLTAGE V $ 

CHARGING MO 0(1 

I w -A| m . CURRENT 
FUNCTIONAL FORMS 

DARK CONDITION • NO PHOTOEMISSION 


SYMBOL A FORMULA 

1/2 


CURRENT DENSITY DUE TO 


ELECTRON BEAM 


SECONDARY ELECTRONS 


BACKSCAmRED ELECTRONS 


LEAKAGE THROUGH INSULATOR 


SURFACE CHARGING 


Figure I - I D model ground test environment 


FLOATING ALUMINUM PLATE 


3 6 i 

s 



0 2 4 6 B 10 

TIME, stc 

Figure 2 - Comparison with piper i mental data 


mv’tonUCIBlLlI Y OF THK 
, PAGE IS POOF 


t 


current, «a 


• CURRENT SllVf RfD TlFlOA 

• VOlTACf r DATA 



hqurtf 3 • Comparison with eiperiroental data 5 mil silvered Teflon 
samples, 100 cm^ area 



Figure 4 Companion with eiperimental data Discharc* transients; silvered 
Teflon samples, 16 kV beam 


SURf ACE VOLTAGE. 


-£ *7C- 


fCXilllBRIUM VOLTAGES 
COMPOSITE SAMPLE 



O ALUMINUM DATA 

a SILVERED TEFLON DATA 

SILVERED TEFLON PREDICTION 

ALUMINUM PREDICTION 


COMPOSITE SAMPlf : S kV BEAM 



U» 


Figure 5 ■ Silvered Itllon/aluminum composite sample floating 
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Figure 7. - Provisional specification for satellite time in a geomagnetic sub- 
storm environment (ref, 11). 
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Figure 8. - Predictions of charging silvered Teflon and aluminum samples. 
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